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Abstract 

CN ' We calculate the branching ratios and CP violating asymmetries of the four B -^ Kt]^'' decays 

CL|' in the perturbative QCD (pQCD) factorization approach. Besides the full leading order contri- 

butions, the partial next-to-leading order (NLO) contributions from the QCD vertex corrections, 
the quark loops, and the chromo-magnetic penguins are also taken into account. The NLO pQCD 

fS| I predictions for the CP-averaged branching ratios are Br{B~^ -^ K~^r]) ~ 3.2 x 10~^, Br{B — > 

K^r]') « 51.0 X 10-6, Sr(S° -^ K^rj) « 2.1 x lO"*^, and Br{B° -^ K^rj') « 50.3 x 10-^. The 

,^ ■ NLO contributions can provide a 70% enhancement to the LO Br{B -^ Kr]'), but a 30% reduc- 

tion to the LO Br{B -^ Kf])^ which play the key role in understanding the observed pattern 
^^, of branching ratios. The NLO pQCD predictions for the CP-violating asymmetries, such as 

^ ; AcpiK^r]') ~ 2.3% and AcpiK^r]') ~ 63%, agree very weh with currently available data. This 

means that the deviation A5 = A^^p{K^r]') — sin 2/3 in pQCD approach is also very small. 
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I. INTRODUCTION 



The B -^ Kr]^'^ decays are very interesting two-body charmless hadronic B meson 
decays. In 1997, CLEO collaboration firstly reported unexpectedly large branching ratios 
for B -^ Kr]' decays pj. Eleven years later, three of the four B — > Kt]^'^ decays have been 
measured with high precision. The world averages as given by HFAG [2| are the following 
fin unit of 10"^) 



Br{B^ - 


-^K^r]) = 2.7 ±0.3, 


Br[B^ - 


, K^n') = 70.2 ±2.5 


Br{B° - 


-^K%) < 1.9, 


Br{B^ - 


-^ K^r]') = 64.9 ±3.1 



A'SUB^ - 


-. K^T]) -- 


= -0.27 ±0.09, 


A'jp{B^ - 


. K^r^') - 


= 0.016 ±0.019, 


A'^,'p{B^ - 


-. K^i) - 


= 0.09 ±0.06, 


A'S'piB^ - 


-. K'v') - 


= 0.61 ±0.07, 



(1) 

From above data one can see that: (a) the measured Br{B -^ Kr]') are much larger 
than the early standard model (SM) expectations, i.e., the so-called fc?7'-puzzle; and (b) 
the large disparity between the branching ratios for B -^ Krj' and B —>■ Krj decays: 
Br{B -^ Kt]') > Br{B -^ Kr]). 

Besides the branching ratios, the CP violating asymmetries for B^ -^ K^^i]^'^ and 
B^ -^ K^i]^'^ decays have been measured very recently [2, l3|: 



(2) 

(3) 

It may be noted that the average of the measured A^^p{B^ —>■ K^t]') is now more than 8cr 
away from zero, so that CP violation in this decay is well established; while Aqp{B^ -^ 
K^r]') is not conflict with zero as expected in the SM. The data for Aqp{B'^ -^ K^r]^'^) 
have less precision, but are consistent with general expectations. 

The measurements of time-dependent CP asymmetries in B^ meson decays, such as 
B^ — > J /^>K^ via h —>■ ccs "tree" transition and B^ — > K^rj' via b —>■ sqq penguin 
transition, have provided crucial tests of the mechanism of CP violation in the SM. Within 
the SM the mixing induced CP violating asymmetry Aqp{B^ — > K^rj') = —rjfSf should 
be comparable with sin 2(3 = 0.685 obtained from the tree dominated B^ -^ J/'^K^ decay, 
this point has been confirmed by the data in Eq. ([3]). 

In the SM the decay B — > Krf^'^ is believed to proceed dominantly through gluonic 
penguin processes!^, [sf and has been evaluated by employing various methods [g, |7|, |8|, 
y, [lO, lul, ll2|, [l3|, [ij, [l5|, [l6|]. Although great progress have been made during the past 
decade, but the predictions for Br{B -^ Kr]') from both the QCD factorization (QCDF) 



approach [ij, [ITJ and the perturbative QCD (pQCD) approach |16l . Il8l | in the Feldmann 



20 are smaller than the data. 



Kroll-Stech (FKS) mixing scheme oi rj — rj' system |19l . 

For the pattern of branching ratios in Eq. ([1]), many possible solutions have been 
proposed. These include, for example, 

(a) Conventional b -^ sqq with constructive (destructive) interference between the 
uu, dd and ss components of 77' (rj) 



(b) Large intrinsic charm content of rl through the chain h — > sec — >■ 577' \P\ or 
through h —>■ sec — > sg*g* — > s{ri,7]') due to the QCD anomaly [Sj]; 

(c) The spectator hard-scattering mechanism through the anomalous coupling of 

gg^v'uMM; 



(d) A significant flavor-singlet contribution [lO|, [l4 

(e) A strong penguin b ^ sg enhanced by new physics [l2|, [l3 . 

But the data of branching ratio in Eq. ([T]) are still not completely understood. For the 
CP violation of -B — *> Kr]^'^ decays, the theoretical studies is still under way. 

In Ref. [l6|, the authors calculated the branching ratios of B ^ Kr]^'^ decays by 
employing the pQCD approach at leading order. They considered the large corrections 
from SU{3) flavor symmetry breaking as well as the possible gluonic component of rj' 
meson, but their prediction for Br{B^ -^ K^rj') ( Br{B^ — > K^rj) ) is much smaller ( 
larger) than the measured value. 

A sizable gluonic content in rj' meson may provide a large enhancement to the decay 



rate of i? — > Kr]'. In Ref. [2l[, the authors examined the possible gluonic contribution to 
the i? — »• 77' transition form factor and found that such contribution is constructive with 
those from quark-content of rj', but numerically very small and can be neglected safely. 
This point has also been confirmed by the QCD sum-rule analysis |22l| 

In the quark-fiavor mixing scheme, the physical rj and rj' meson are linear combinations 
of flavor state ?7g = {uu+dd)/y/2 and r]s = ss with the "mass" of irigq and m^s respectively. 
In Ref. [23|], the effect of a large chiral scale ml = m'^g/{2mq) with q = {u,d) for the 
meson rjq has been evaluated although we do not know which mechanism is responsible to 
achieve a large value of m^g. When one uses m^g = 0.22 GeV [23| instead of its generally 
accepted value of rrigq = 0.11 GeV, a larger B -^ Krjg decay amplitude can be obtained. 
Consequently, the LO pQCD predictions for Br{B -^ Krj') become consistent with the 
data. 

In Ref. [2J], the authors examined the possible way to increase the value of mqq. 
They found that few-percent violation of Okubo-Zweig-Iizuka (OZI) rule can enhance 
rriqq few times, which then leads to the consistency of the LO predictions with the data 
for B — >■ Kri'^''> decays. 

Besides the possible mechanisms mentioned above, we here consider a new and nat- 
ural solution: the effects of the next-to-leading order (NLO) contributions in the pQCD 



approach. As shown in Ref. [25|, the NLO contributions to i? — > Kti decays can play the 
key rule to explain the so-called "i^vr"— puzzle. We expect here the NLO contributions 
could help us to resolve the ^^Kt]'" —puzzle. 

For the CP asymmetries of B^ —>■ K^rj', the deviation ASf = —i]fSf — sin 2(3 has been 
estimated, for example, in the QCDF approach [l5|, |26| and the soft coUinear effective 
theoryJ27l|. The resultant bound is \ASf\ < 0.05. Since the source of the CP violation in 
the pQCD approach is very different from those in the QCDF/SCET approach, we here try 
to calculate the CP asymmetries oi B -^ Kt]^'^ decays by employing the pQCD approach 
at LO and NLO level, to check if we can accommodate the data of CP asymmetries. 

In this paper we will calculate the next-to-leading order contributions to the branching 
ratios and CP violating asymmetries of the four B — > Kr]^'^ decays. We firstly calculate 
the decay amplitudes of the B — > Kr]^'^ decays by employing the pQCD factorization 



approach at the leading order (LO), as have been done in previous studies for other 



two-body charmless B meson decays [28|, |29|, |30|, l311- And then we evaluate the NLO 
contributions to these decays. 

The NLO contributions considered here include: QCD vertex corrections, the quark- 
loops and the chromo-magnetic penguins. We wish that they are the major part of the 
full NLO contributions in pQCD approach [25|. Of course, remaining NLO contributions 
in pQCD approach, such as those from factorizable emission diagrams, hard-spectator 
and annihilation diagrams, should be calculated as soon as possible. 

This paper is organized as follows. In Sec. II, we give a brief review about the pQCD fac- 
torization approach. In Sec. Ill, we calculate analytically the relevant Feynman diagrams 
and present the various decay amplitudes for the studied decay modes in leading-order. 
In Sec. IV, the NLO contributions from the vertex corrections, the quark loops and the 
chromo-magnetic penguin amplitudes are evaluated. We calculate and show the pQCD 
predictions for the branching ratios and CP violating asymmetries oi B ^ Kr]^'^ decays 
in Sec. V. The summary and some discussions are included in the final section. 



II. THEORETICAL FRAMEWORK 

A. Theoretical framework 

In the pQCD approach, the decay amplitude is separated into soft ($AfJ, hard 
( H{ki,t) ), and harder( C{Miy) ) dynamics characterized by different energy scales 
(Aqc'j),^, m^, M14/) [l8|. The decay amplitude A{B — >■ M2M3) can be written conceptu- 
ally as the convolution, 

A{B ^ M2Ms) --^ fd%d%d% TT[CmB{ki)<^M,{k2)^M^{h)H{h,k2,h,t)], (4) 

where /cj's are momenta of light quarks included in each meson, and Tr denotes the trace 
over Dirac and color indices. C{t) is the Wilson coefficient evaluated at scale t. In 
the above convolution, the Wilson coefficient C{t) includes the harder dynamics at scale 
higher than Mb and describes the evolution of local 4-Fermi operators from mw ( the 
W boson mass) down to t ~ 0{^/AMb) scale, where A = Mb — rrib. The function 
H{ki, /c2, ^3, t) describes the four quark operator and the spectator quark connected by a 
hard gluon whose g^ is in the order of AMb, and includes the (9( y AM^) hard dynamics. 
Therefore, this hard kernel H can be perturbatively calculated. The function $A/i is the 
wave function which describes hadronization of the quark and anti-quark in the meson 
Mi. While the hard kernel H depends on the processes considered, the wave function $Mi 
is independent of the specific processes. Using the wave functions determined from other 
well measured processes, one can make quantitative predictions here. 

Since the b quark inside the B meson is rather heavy, we consider the B meson at rest 
for simplicity. It is then convenient to use light-cone coordinate {p'^,p~ , prp) to describe the 
meson's momenta: p^ = 4f (p'^ ±p^) and px = {p^jP"^) ■ Using the light-cone coordinates 
the B meson momentum Pb and the two final state meson's momenta P2 and P3 (for M2 
and M3 respectively) can be written as 

Pb = -^{1,1,0t), P2 = -^{l~rlrlOT), P3 = ^(r3^ 1 - 4 Ot), (5) 



where r^ = rrii/MB- m2 and m^ are the mass of the two final state mesons. For the case 
oi B —> PP decays, r2 and r^ are small and could be neglected safely. 

Putting the anti-quark momenta in B, M2 and M3 meson as ki, /c2, and k^, respectively, 
we can choose 

fci = (xiP+,0,kiT), k2 = {x2P^,0,k2T), ks = {0,XsPi,ksT). (6) 

Then, the integration over k^ , k2 , and k^ in eq.(jl]) will lead to 



A{B -^ PV) ~ dxidx2dxsbidbib2db2bsdbs 

■Tr [CmBixub^)^MAx2,b2)<^M,{^3,bs)H{x^,h,t)St{x,)e-''^'^] , (7) 

where bi is the conjugate space coordinate of kir- The large logarithms (\iamw/t) coming 
from QCD radiative corrections to four quark operators are included in the Wilson coeffi- 
cients C(t). The large double logarithms (In Xi) on the longitudinal direction are summed 
by the threshold resummation, and they lead to St{xi) which smears the end-point singu- 
larities on Xi. The last term, e~^^^\ is the Sudakov form factor which suppresses the soft 
dynamics effectively |l8| . 



B. Effective Hamiltonian and Wilson coefficients 

For the studied B — > KT]^decajs, the weak effective Hamiltonian H^ff for b 
transition can be written as 3^ 

* g=u,c (, i=3 

where Gp = 1.16639 x lO^^GeV'"^ is the Fermi constant, and Vij is the CKM matrix 
element, Ci{fi) are the Wilson coefficients evaluated at the renormalization scale n and 
Oi{fi) are the four-fermion operators. For the case of 6 — *■ (i transition, simply makes a 
replacement of s by c? in Eq. ([8]) and in the expressions of Oi{^) operators, which can be 
found easily for example in Refs. 30|, ImI, l32 • 



In PQCD approach, the energy scale "t" is chosen as the largest energy scale in the 
hard kernel H{xi, bi, t) of a given Feynman diagram, in order to suppress the higher order 
corrections and improve the reliability of the perturbative calculation. Here, the scale "t" 
may be larger or smaller than the rrij, scale. In the range of t < mi, or t > rrih, the number 
of active quarks is Nj = 4 or Nf = 5, respectively. For the Wilson coefficients Ci{fi) and 
their renormalization group (RG) running, they are known at NLO level currently |32| . 
The explicit expressions of the LO and NLO Ci{my/) can be found easily, for example, in 
Refs. fiii. 



When the pQCD approach at leading-order are employed, the leading order Wilson 
coefficients Ci{mw)i the leading order RG evolution matrix f/(t, m)^'^'> from the high scale 
m down to t < ??i ( for details see Eq. (3.94) in Ref. [3^), and the leading order as{t) are 
used: 

47r 

«s(t) = ^ , r,^,.2 T, (9) 



where (3o = (33 - 2Nf)/3, AgJ.^ = 0.225G'e\/ and A^J.^ = 0.287 GeV. 

When the NLO contributions are taken into account, however, the NLO Wilson coef- 
ficients Ci{mw), the NLO RG evolution matrix U(t,m,a) ( for details see Eq. (7.22) in 



Ref. [32|) and the as{t) at two- loop level are used: 



where po = (33 - 2A^j)/3, /3i = (306 - 38iV/)/3, AgJ,^ = 0.225 GeV and AgJ,^ = 0.326 
GeV. 

From the general knowledge, the hard scale t must be much larger than Aqcd ~ 
0.2 GeV in order to guarantee the reliability of perturbative calculations. In previous 
calculations based on the pQCD approach /io = 0.5 GeV is chosen as the lower cut-off of 
the scale t. In our opinion, it is indeed too low, because it may be conceptually incorrect 
to evaluate the Wilson coefficients at scales down to 0.5 GeV [33]. The explicit numerical 



checks as done in Ref. [3J] also show that (a) the Wilson coefficient Ci(0.5) is close to 
—1 and clearly too large in size! (b) the values of the Wilson coefficients C3,4,5,6(/^) at 
/i = 0.5 GeV are about four to seven times larger than those at /i = 1.0 GeV; and (c) the 
/iQ— dependence of all Wilson coefficients become relatively weak for /iq > 1.0 GeV. We 
therefore believe that it is reasonable to choose /ig = 1.0 GeV as the lower cut-off of the 
hard scale t, which is also close to the hard-coUinear scale yArriB ~ 1.3 GeV in SCET. 
In the numerical integrations we will fix the values Ci{t) at Ci(l.O) whenever the scale t 



runs below the scale /xq = 1.0 GeV [SJ, [35 



C. Wave functions 

Since the b-quark is much heavier than the up or down quark, the B meson is treated 
as a very good heavy-light system. Although there are in general two Lorentz structures in 
the B meson distribution amplitudes, they obey to the following normalization conditions 

'^''M^.)-^^, /^MkO-o. (11) 



(2^)4--v -^ 2v/2iV;' J (27r) 

However, it can be argued that the contribution of 4>b is numerically small [36], thus its 
contribution can be numerically neglected. In this approximation, we only consider the 
contribution of Lorentz structure 

$B = -^i^B + msHM^i), (12) 



with 



(l)B{x,b) = Nbx'^{1 — xYexp 



Ml x^ 1, ,,2' 



2uj^ 2 



(13) 



where u)b is a free parameter and we take Ub = 0.4 ± 0.04 GeV in numerical calculations, 
and Nb = 101.445 is the normalization factor for ut, = 0.4. 

6 



The Kaon mesons are treated as a light-light system. The wave function of K meson 
is defined as 37 



$x(P,x,C) 



1 



V2iv; 



=75[f'0^(x)+m^0^(x) + Cm, 



V ■ n] 



'K 



X 



(14) 



c 



where P and x are the momentum and the momentum fraction of K, respectively. The 
parameter C, is either +1 or —1 depending on the assignment of the momentum fraction 

X. 

For ri^''> meson, the wave function for rfq components of r]' meson are given as 



%iP,xX) = 



ym 



=75 



c 



x) + ml(p^{x) + C"^o(# - ^ ■ n)<plix) 



(15) 



where P and x are the momentum and the momentum fraction of r^g, respectively. We 
assumed here that the wave function of rjq is same as the vr wave function. The parameter 
C, is either +1 or —1 depending on the assignment of the momentum fraction x. The 
rjs = ss component of the wave function can be defined in the same way. 

The expressions of the relevant distribution amplitudes (DA's) of K meson are the 
following [37|: 



<f>i{x) 

0^(X) 

0^(x) 



/. 



K 



Qx(l — x) 



Kn'^11 



Kr^'i/'i, 



Kri'i/^i 



1 + afCl'\t) + 4C^'\t) + 4Cl'\t) 






9 

^3^3 + — P^(l + 6af) 



/ 



K 



1 + 6(5r73 - \riz^z - ^pl " ^P^«f )(1 " 10:^ + lOx^) 



(16) 

Cl'\t) \ (17) 

18) 



2V2iVc 
with the mass ratio pk = fnK/rrioK- The Gegenbauer moments can be given as |37] : 



K 



.K 



K 



0.2, a^ = 0.25, a2 = -0.015. 



(19) 



The values of other parameters are 773 = 0.015 and uj = —3.0. At last the Gegenbauer 
polynomials C^{t) are given as: 

Cl^\t) = li^t' - 1), Cy\t) = ^(3 - sot' + 35t^), 



y3/2 



Ct^\t) 



■y3/2 



3t, 

15 



Cf(t) = |(5t^-1), 



C2^\t) = — (l-14t^ + 21t^ 



(20) 



with t = 2a; — 1. 

In the quark-flavor mixing scheme, the physical states r] and r]' are related to the flavor 
states rjg = {uu + dd)/y2 and rjs = ss through a single mixing angle 0, 



cos I 
sin< 



— smi 
cosd) 



Vs 



Fi{(f)){uu + dd) + ^2(0) ss 
F[{<p)luu + dd) + F!2^{<p) ss 



(21) 



with 



Fm 



coscp 
sin(/) 



V2 



F2{(f)) = -siiii 
F^{(j)) = cos(f). 



:22i 



as 



The relation between the decay constants (/^, /*, /^,, /^,) and (/,, /«, ) can be written 

(23) 



cos (p — sm ( 
sin 6 cos (h 



U 

/. 



fq fs 
J rj J Tj 
f? f s 

The chiral enhancement mg and rriQ associated with the two-parton twist-3 r]q and t], 
meson distribution amphtudes have been defined as [2 

i2 1 



iTin = T-^ = - — Imt cos'^ (f) + mi, sin^ (j) ^fTTi^/ — m^) cos(/)sin(/)l, (24) 



■rrir, 



2m„ 



mz 



2m^ 



r 2 2 J I 2-2 

- — m„ cos + m„, sm . 
2m„ ' ' 



V^/s/ 2 
Jq 



r 2 2 J I 2-2 

m / COS + m sm , 

2mo ' ' 



'9 

V2/, 



/s 



9 / 2 



)cos0sin0], (25) 



m^ ) cos sm 



by assuming the exact isospin symmetry ruq = mu = rrid. The three input parameters 
fq, fs and (j) have been extracted from the data of the relevant exclusive processes [l9|: 

/, = (1.07 ± 0.02)/^, /, = (1.34 ± 0.06)A, = 39.3° ± 1.0°, (26) 

The distribution amplitudes 0J^'^'^ represent the axial vector, pseudoscalar and tensor 
component of the wave function respectively [37| . They are given as: 



<(^) 



C(^) 



2v/2iV: 
+al''Cl^\2x - 1 



riqr^?./2 



Qx{l -x) 1 + atCl^\2x - 1) + al'C'^'\2x - 1 



»?9^3/2. 



(27) 



2^/W^ _ 



l + (30r/3--pt)C2' (2x-l 



1/2. 



2^^g 



9 



2 

20 ^'^^^ 



-3<'r73W3 + — p^ fl + 6a^'; 



C(^) 



/^ 



:(1-2X) 



2v/2Ar: 

■(l-lOx+lOa;^)] 



6 Sr/s 



'7? 



0, a^' 



Cl/'(2x-l) 


5 


(28) 


1 7 2 
-2^3^3 20^^^. 


3 2 .,\ 
- 5P..-2' J 




(29) 


0.44 ± 0.22, a 


^9 

4 


= a\ = 0.25, and the 



where p^^ = 2mq/mqq,a{ 
Gegenbauer polynomials C'^{t) have been given in Eq. (1201) . As to the wave function and 
the corresponding DA's of the ss components, we also use the same form as qq but with 
some parameters changed: p^^ = 2ms/mss, a^" = al'' for i = 1,2,4. 

The transverse momentum fc^ is usually converted to the b parameter by Fourier 
transformation. The initial conditions of leading twist (f)i{x), i = B,K,ri,ri', are of non- 
perturbative origin, satisfying the normalization 



diix, b = 0)dx = — ;= 
'0 2VQ 

with fi the meson decay constant. 



fi 



(30) 



8 



B 




(a) 



K\ri^'^] B 



(b) 




(c) 



(d) 






ri<^'\K\ 



oc^ 




(g) (h) 

FIG. 1: Feynman diagrams which may contribute to the B 



Krj^'' decays at leading order. 



III. DECAY AMPLITUDES AT LEADING ORDER 



In the pQCD approach, the Feynman diagrams as shown in Fig. [T]may contribute to 
B -^ Kt]^'^ decays at leading order. As mentioned previously, B^ -^ K^r]^'^ decays have 
been studied in Ref. [l6j by employing the LO pQCD approach. In this section, we firstly 
calculate the LO decay amplitudes for four B -^ Ki]^'^ decays, but in a rather different 
way to treat the Feynman diagrams from that in Ref. [16| . 

At the leading order in pQCD approach, there are three type diagrams contributing to 
the B -^ Kt]^'^ decays, the factorizable emission diagrams, the hard-spectator diagrams 
and the annihilation diagrams, as illustrated in FigJH From the factorizable emission 
diagrams, the corresponding form factors can be extracted by perturbative calculation. 
First, we consider the B -^ Krj decay modes, and then extend the calculation to B ^ Krj' 
decays. 

For the usual factorizable emission diagrams 1(a) and 1(b) with the B -^ K transition, 
i.e., it is the K meson pick up the spectator quark, the operators Oi, O2, O34 and Og^io 



are {V — A){V — A) currents, the sum of the individual amphtudes is given as 

FeK = -^'nGFCFm^B I dxidx2 I bidbib2db2(l)Bixi,bi) 
v2 Jo Jo 

X { [(1 + X2)0:^(X2) + (1 - 2X2)rKi(pKix2) " 0x(^2))] " ^e(ta)/ie(Xl, X2, 61, ^2) 

+2rK(^^{x2) ■ E,{QK{x2, xi, 62, &i)} , (31) 

where tk = m^ /mB with itiq is the chiral scale; Cp = 4/3 is a color factor, and X2 = 
1 — X2. The evolution function Ee{t) and hard function h^ are displayed in Appendix 
lAl In the above equation, we do not include the Wilson coefficients of the corresponding 
operators, which are process dependent. They will be shown later in the expressions of 
total decay amplitude. 

Also for diagrams 1(a) and 1(b), the operators O^^ and O^^s have a structure of (V — 
A){V + A) currents. In some decay channels, some of these operators contribute to 
the decay amplitude in a factorizable way. Since only the axial-vector part of {V + A) 
current contribute to the pseudo-sealer meson production, {K\V — A\B){r]^'^\V + A\0) = 
-{K\V - A\B){r]^'^\V - A\0), that is 

Ff^ = -F,K . (32) 

In some other cases, we need to do Fierz transformation for those operators to get right 
color structure for factorization to work. In this case, we get {S — P){S + P) operators 
from {V — A){V + A) ones. For these {S — P){S + P) operators, the corresponding decay 
amplitude is 

1 f^ r^ /'oo 

F^3 = —=nGFCF'm'i, / dxidx2 / bidbib2db2<i>B{xi] 
V2 Jo Jo 

X {"^V [0k(^2) + ^i^((2 + X2)<Pk{^2) + a;2</>^(x2))] ■ Ee{ta)he{Xi, X2, 61, 62) 

+2rKrr,(pK{^2) ■ Ee{Qhe{x2,Xi,b2,bi)} . (33) 

where r^ = mQ/rriB, and ml = rnl'' is the chiral scale defined in Eq. (12^ . 

For the non- factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 62 can be performed using 5 function 5{b^ — 62), leaving only 
integration of 61 and 63. For the {V — A)(y — A) operators, the result is 

16 [^ f°° 

MeK = -^T^GpCFTn^B / dxidx2dx3 / bidbibsdbs(j)B{xi,bi)(f)^{x3) 
V 3 Jo Jo 

X { [-rKX2 (0^(X2) + 0^(X2)) + (1 - X3)(j)i{x2)] 

■E'^{tb)hn{Xi, X2, 1 - X3, bi, 63) 
+ [-(a;2 + 3:3)0^(^2) + rKX2 (</>^(^2) - </>x(^2))] 

■K{t'h)hn{.Xi, X2, X3, 61, 63)} , (34) 

where 0,, denotes 0^^ or 0^^. 

There are two kinds of contributions from (V^ — A){y + A) operators: M^ and M^, 
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corresponding to the (V — A){V + A) and {S ~ P){S + P) type operators respectively: 

v3 Jo io 

■{[(1-X3)0^(X2) (0^(0:3) -0J(X3)) 

+r;^(l - 3:3) (0^(X2) + 0^(0:2)) (0j'(5;3) - 0^^(2:3)) 

+rKX2 (0^(X2) - 4>k{x2)) (0j(^3) + ^^(^s))] 
■E'^{tb)hniXi, X2, 1 - ^3, 61, 63) 

- [ 3:30:^(^2) (0j(^3) + ^Jfe)) 

+r2X3 (0^(X2) + 0^(2:2)) {(P^ixs) + 05^(0:3)) 

+^23^2 (0:^(^2) - ^^fe)) (0jfe) - 0j(^3)) ] 
■Ki't'bJhniXi, X2, X3, 61, 63)} , (35) 

16 [^ f°° 

M[j^ = —=7rGFCFm'^ / dxidx2dx3, / 61^6163^63 0B(a;i, 61)0^^(3:3) 

v3 io io 

• { [-(1 + X2- ^3)0^(^2) + a;2r7^ (^^(^2) - 0]^(^2))] 

■E'^(tb)hn{Xi, X2, 1 - X3, 61, 63) 
+ [X30^(X2) -X2rK (0^(^2) +^^(^2))] " ^e(t'J/l„(Xi, X2, X3, 61, 63) } . (36) 

For the non-factorizable annihilation diagrams 1(e) and 1(f), again all three wave 
functions are involved. Here we have two kinds of contributions: M^ = 0, MaK and 
M^ describe the contributions from the {V — A)(y — A) and {V — A)(y + A) type 
operators, respectively, 

1 fi /*! /»oo 

MaK = —^TrGFCF'm'^ / dXidX2dx3, / 61^6163^63 0^(3^1, 61) 

v3 Jo Jo 

■ { [(1 - X2)0^(x2)0;^(a;3) + VKT^il - X2) (0^(X2) + 0^(a;2)) (0^(0:3) - 0^(2:3)) 

+rKrr,X3 (0x(^2) - 0x(^2)) (0j'(S3) + ^^(^s))] " E'^{tc)hna{Xl, X2, X3, 61, 63) 
- [0:30^(^2)0;? (^3) + 4ri^r^0^(x2)0j'(S3) - rxrr,(l - Xg) (0x(^2) + 0^(^2)) 

■ (0j'(^3) - (P^ixs)) - rKrr,X2 (0x(^2) - 0x(^2)) (0j'(%) + 05^(^3))] 

X ■E^(t'J/i„a(xi,a;2,X3,6i,63)} , (37) 

^cS = ^pTrG'irCiT'm^ / dxidx2dx3, / 61^6163^63 0B(a;i, 61) 
v3 Jo Jo 

X { [-(1 - X2)rK(t)^{x3) (0^(^2) + 0^(^2)) + r^a;30^(x2) (0J(X3) - 0j'(^3))] 

x£'^(tc)/i„a(a:i,a:2,X3,6i,63) - [(^2 + l)r/^0;?(x3) (0]^(x2) + 0]^(S2)) 

+r,(x3 - 2)0^(^2) (0j(x3) - 0^(2:3))] ^:(t'J/i„„(xi,a;2,X3,6i,63)} . (38) 

The factorizable annihilation diagrams 1(g) and 1(h) involve only K and r]'^''^ wave 
functions. There are also three kinds of decay amplitudes for these two diagrams. FaK-, 
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i^fi and F^l 



Q /•! poo 

V^ ^0 JO 



■Ea{td)ha{Xs,l - X2,b3,b2) 
+ [X30^(X2)0;^(X3) + 2r^rK0^(x2) (0J'(X3) + 0j'(^3)) 
+2r^r/^X30]^(x2) (0J'(X3) -0J'(X3))] ■Ea(t'^)/la(l-X2,X3,62,&3)} , (39) 



.P2 16 



F„^ = —=7iGFCpm^ / dx2dx3 / b2db2hdb3 
v2 Jo io 

■ { [r^(l - X2) (0^(X2) - 0^(X2)) <(X3) + 2r,0^(x2)0j'(x3)] 

■Ea{td)ha{x3, 1 - X2, &3, ^2) 

+ [2r^0^(x2)0;^(x3) +X3r^0^(x2)(0j'(x3) +0^(0:3))] 

■Ea(t'J/l,(l-X2,X3,&2,&3)} • (40) 

The evolution function Ei(tj) and hard function hi appeared in Eqs. (13311401) are given 
exphcitly in Appendix El 

If we exchange the K and i]^'^ in Fig. 1, the corresponding decay amphtudes for new 
diagrams will be similar with those as given in Eqs. fl31H40l) . since the K and t]^'-* are all 
pseudoscalar mesons and have the similar wave functions. The decay amplitudes for new 
diagrams, say Fe^, Ff^^'^"^, M^r,, M^^'^"^, Mar^.M^^, Fa^,F^^^'^^, can be obtained from 
those those as given in Eqs.( l3Tll40l) by the following replacements 



K ^ (P^C) . (pK ^ (P^C) . (pK ^ (P^C) . ^K ^ r^c) . (41) 

For B^ -^ K^T] decay, by combining the contributions from all possible configuration 
of Feynman diagrams, one finds the total decay amplitude with the inclusion of the 
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corresponding Wilson coefficients as follows 
M{K^7]) = < K^r]\H^ff\B^ >=F, 



eK 



1 



^«ct2 - 6 ( 2a3 - 2a5 - -ay + -ag 



fq 



-6 ( as + 04 - ^5 + 2*^7 - 2^9 - ^ctio ) /^ 



1 

— ( 

2 



- [FakF2i<p) + Fa^F^m 6 U4 - \aio 






^6 — 2*^8 ) /s 



FM 



-it (C, + C4 - iCg - iCio") i^2(<^)} 

-Mi|6 (^2C6 + ic8^Fi(</)) + (C6-^C8)F2(</)) 

where ^^ = Vu*ftVus, 6 = ^fe^s, and -Fi(0), -^2(0) are the mixing factors as given in Eq. ( l22l) . 
The coefficients Oj in Eq. fH2]) are the combinations of the Wilson coefficients Cj, and 
have been defined as usual 



(42) 



di = C2 + — , aa = Ci + — , 

ai = Ci-\ ^, for i = 3,5, 7, 9, 

at = Ci + ^, for z = 4,6,8,10. 



(43) 
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Similarly, the decay amplitude for B^ -^ K^rj can be written as 



MiK+ri) = < K+ri\Heff\B° > 



eK 



iua2 - it 2a3 - 2a 



1 1 



fq 



-6 ( % + 04 - 05 + -Or - -flg - -aio ) /^ 



- [Fe^Fi(0)/x + (F,,Fi(0) + F„;^F2(0)) /b] 6 (^4 + «io) 



.PP2 



ei^/)^6 ( "6 - -t 



Mfi6 



a 



> 



+M, 



eii' 



'C«C'2 



.^i ( 2C4 + -Cio 



i^l(0) 



-6 C3 + C4 



iCg - iCio ) F,{<P) 



+ [MaKF^icj)) + (Me, + M,,) Fi(0)] [e„Ci - 6 (C3 + Cg)] 

- [Mf^F2(0) + (Mf^i + Mf^i) Fi(0)] 6(C5 + C77) 



P2( 



2C6 + -C8 Fi(0) 



-Mf^^e 



The total decay amplitude for B^ -^ K'^rj' and -B^ 
from Eqs. fl42p and fl44p by the following replacements 



Ce — ■^C's 



i^2(0) 



(44) 



K^f]' can be obtained easily 



Fi{<P) 



Jr)'i J ri Jr)'i 

F[{<P), F2{(t>) ^ F^{<t>). 



(45) 



IV. NLO CONTRIBUTIONS IN PQCD APPROACH 



A. General discussion 



The power counting in the pQCD factorization approach |25j is different from that 
in the QCD factorization IJ, |r7|. When compared with the previous LO calculations in 
pQCD [l8|, |30|, l3l| , the following NLO contributions should be considered: 



1. The LO Wilson coefficients Ci{mw) will be replaced by those at NLO level in NDR 
scheme 32|, and the NLO RG evolution matrix U{t,m,a) instead of U{mi,m2)^^\ 
as defined in Ref. 32 , will be used here: 



a 



f/(mi, m2, a) = f/(mi, 7712) + -— -R(mi, 777-2) 

Air 



(46) 



where the function U{mi, 7712) and R{mi, 7712) represent the QCD and QED evolution 
and have been defined in Eq. (6.24) and (7.22) in Ref. 32]. We also introduce a 
cut-off /io = 1.0 GeV for the QCD running of Ci{t) in the final integration. 
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(a) (b) (c) (d) 

FIG. 2: NLO vertex corrections to the factorizable amplitudes. 

2. The strong couphng constant as{t) at two- loop level as given in Eq. (ITOl) will be 
used. 

3. Besides the LO hard kernel H^'^^ag), the NLO hard kernel H^^^a'j.) should be 
included. All the Feynman diagrams, which lead to the decay amplitudes propor- 
tional to a^(t), should be considered. Such Feynman diagrams can be grouped into 
following classes: 

1: The vertex corrections, as illustrated in Fig. [2l the same set as that studied 
in the QCDF approach. 

II: The NLO contributions from quark-loops, as illustrated in Fig. [3l 

III: The NLO contributions from chromo-magnetic penguins, i.e. the operator 
Osg, as illustrated in Fig. HI There are totally nine relevant Feynman diagrams 
as given in Ref . |38|] , if the Feynman diagrams involving three-gluon vertex are 
also included. We here show the first two only, and they provide the dominant 
NLO contributions, according to Ref. |38|]. 

IV: The NLO contributions to the Feynman diagrams (la, lb) corresponding 
to the extraction of from factors, as illustrated in Fig. [51 There are totally 13 
relevant Feynman diagrams, we here show four of them only. 

V: The NLO contributions to the hard-spectator Feynman diagrams (lc,ld), 
as illustrated in Fig. [61 There are totally 56 relevant Feynman diagrams, we 
here show four only. 

VI: The NLO contributions to the annihilation Feynman diagrams (le,lh), as 
illustrated in Fig. [3 we here show only four such diagrams. 

For the last four classes (III- VI), the Feynman diagrams involving three-gluon vertex 
should be included. At present, the calculations for the vertex corrections, the quark-loops 
and chromo-magnetic penguins have been available and will be considered here. For the 
Feynman diagrams as shown in Figs. 5-7, however, the analytical calculations have not 
been completed yet. What we can do here is to include the NLO contributions to the 
hard kernel H. 

B. Vertex corrections 

The vertex corrections to the factorizable emission diagrams, as illustrated by Fig. [2l 
have been calculated years ago in the QCD factorization appeoach|14j. Il5l . llTI . 
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(a) 



(b) 



FIG. 3: Quark-loop amplitudes. 

?/')(A') 








(a) 



(b) 



FIG. 4: Chroino- magnetic penguin amplitudes (Ogg)- There are nine relevant Feynman diagrams 
as shown in Ref. [33]. Here we show the first two only, which provide dominant contribution of 
such diagrams. 



For the emission diagram, there are 4 kinds of single gluon exchange responsible for the 
effective vertex as labeled in Figj2l The contributions from the soft gluons and coUinear 
gluons are power suppressed, that is to say the total contributions of these four figures are 
infrared finite. For charmless B meson decays, these corrections can be calculated without 
considering the transverse momentum effects of the quark at the end-point in collinear 
factorization theorem. Therefore, there is no need to employ the fcr factorization theorem. 
In fact, the difference of the calculations induced by considering or not considering the 
parton transverse momentum is rather small [25], say less than 10%, and therefore can 
be neglected. Consequently, one can use the vertex corrections as given in Ref. 15 



directly. The vertex corrections can then be absorbed into the re-definition of the Wilson 
coefficients ai{fi) by adding a vertex-function Vi{M) to them |15l. Il7 



aj(/x) 



a^{^^) + '^^Cf^V,{M), for ^ = 1, 2; 






47r 

«s(/i) ^ Cj+i(/i) 
47r ^ 3 



V;(M), for J = 3, 5, 7, 9, 



VjiM), for J =4, 6, 8, 10, 



(47) 



B 



\jC- \/{, ^/, ,,\/^ 



y: 



FIG. 5: The four typical Feynman diagrams, which contributes to the form factors at NLO level. 
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B 



X^ 




FIG. 6: The four typical hard-spectator Feynman diagrams, which contributes at NLO level. 




FIG. 7: The four typical annihilation Feynman diagrams, which contributes at NLO level. 

where M is the meson emitted from the weak vertex. When M is a pseudo-scalar meson, 
the vertex functions Vi{M) are given ( in the NDR scheme) in Refs. [l5|, |25|: 



l2ln^~18 + ^J^dx<PUx)g{x), for z = 1-4,9,10, 



V,{M) 



-121n^ + 6 
-6 



2v^ fl 



2V6 fl 



/a 



h^J^'dxKi^)g{l-x) 
Jo dx(f)f^{x)h{x), 



for i 
for i 



5,7, 
6,8, 



(4J 



where /m is the decay constant of the meson M; 4'm{x) and 0^/(a;) are the twist-2 and 
twist-3 distribution amplitude of the meson M, respectively. The hard-scattering functions 
g{x) and h{x) in Eq. fH5]) are: 



9{x) 



h{x) 



2x 



1 



\nx — in 



+ 



21nx 



2Li2{x) -ln^x + 
2Li2{x) - In^ X - (1 + 2iTx) Ins 



-(3 + 2z7r)lnx 
1 - 



[X 



[X 



X) 



(49) 
(50) 



where Li2{x) is the dilogarithm function. As shown in Ref. |25j, the /^-dependence of 
the Wilson coefficients aj(/x) will be improved generally by the inclusion of the vertex 
corrections. 



C. Quark loops 

The contribution from the so-called "quark-loops" is a kind of penguin correction with 
the four quark operators insertion, as illustrated by Fig. [3l In fact this is generally called 
the BSS mechanism |3 91] ■ which provide the strong phase needed to induce the CP violation 
in QCDF approach. We here include quark-loop amplitude from the operators Oi^2 and 
O3-Q only. The quark loops from Oy^io will be neglected due to their smallness. 
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For the 6 — > s transition, the contributions from the various quark loops are given by: 



H. 



eff 



Gf 



as(/i) 



E E ^l^^^^''*^^^'^'"' ^') ^'^^ (^ - ^^) ^"^) (?"7''T»g') , (51) 



q=u,c,t qi 



where /^ is the invariant mass of the gluon, which attaches the quark loops in Figl31 The 
functions C"^(/i, /^) are written as 



c^i^iX' 



G'^ifij')- 



C2(/x), 



(52) 



for q = u,c and 



G^'\fi,f 



C3(/i)+ Yl Giq")if^,f)[C,ifx) + G,ifx)]. (53) 

q//=u,d,s,c 



The function G^^\fj., P) for the loop of the massive q{q = u, d, s, c) quark is given by |25| 



G(«)(/i,/2) = _4 / clxx{l-x)ln^ ^^ ^' 



/^^ 



(54) 



rrir 



is the possible quark mass. The explicit expressions of the function G^'^\n,l'^) after 
the integration can be found, for example, in Ref. [25l |. 

It is straightforward to calculate the decay amphtude for FigjS^ and [3]d. We find two 
kinds of topological decay amplitudes: 



m: 



(g) 



K„ — ■i=Cpm% / dxidx2dx3 / hidhib2dh2(t)B{xi,hi) 

v6 Jo Jo 

■ { [(1 + X2) 0^(a:2)<(x3) + Ti, (1 - 2X2) (0^(X2) - 0^(X2)) <(X3) 

+2r^>;^(x2X(a;3) + 2rKr^, ((2 + X2) ^^(xs) + X20^(x2)) <(a;3)] 

■E'^'i\t^,l^)K{x2,xiMM) 
+ [2rK0x(a;2)0;t(^3) +4r;^r^,0^(x2)0j(x3)] 

■E(«)(t;,0/ie(xi,X2,6i,62)}, (55) 

for B ^ K transition, and 






-=Cp7n'^ / dxidx2dx3 / 61^6162^^2 0_b (2^1 , ^1) 

v6 Jo Jo 

• { [(1 + X2) <(X2)0^(X3) + r, (1 - 2X2) (0J;(X2) - 0J(X2)) ^^(xa) 

+2^/^01 (^2)0x(^3) + 2r^rK ( (2 + Xa) ^J (^2) + Xa^J^fxa) ) ^xfe) 



+ 



•E('')(t„/2)/i^(a;2,xi,62,6i) 
2r,,0j'^(x2)0:^(x3) + 4r^ri^0j'^(x2)0:^(x3) 

■E('')(t;,0/ie(Xl,X2,6i,62)}, 



(56) 
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for S — ^ ?7 transition. Here r^ = ml/mB and r.^^ = mf^/mB. The evolution factors in 
Eqs. dSSEnD take the form of 

EW(t,/2) = C^'i){t,l')al{t)-eM-Sab]. (57) 

with the Sudakov factor Sab and the hard function he{xi,X2,bi,b2) as given in Eq. ( 1A9I) 
and Eq. (IA2p respectively, and finally the hard scales and the gluon invariant masses are 

tg = max{^/a^mB,^/x^mB,^/iT^^X2ja^mB,l/bl,l/b2)^,, 

t'g = max{^/x^mB,^/a^mB,^/\x^^^x^\mB,l/bl,l/b2), (58) 

/2 = (1 _ X2)xsm% - |k2T - kgTp ?^ (1 - X2)x3m%, 

V^ = (x3 - xi)m| - |kiT - ksTp ^ (x3 - xi)m|. (59) 

For i? -^ Kf]' decays, we find the same decay amplitude. Finally, the total "quark- 
loop" contribution to the considered B — ^ Ki]^'^ (K = K'^^K^) decays can be written 
as 

M(,^y = < Kr^\ntf,\B >= ^ 5^ \ [m^^ F2(0) + M^^), Fi(0)] , (60) 

* q=u,c,t 

M^^i^, = < KViO >= ^ $: A, [M^t ^2(0) + M^:!, F{(0)] , (61) 

^ q=u,c,t 

where Xq = VgbV*],. The mixing parameters -Fi(0), F[{(j)), ^2(0) and -^2(0) have been 
defined in Eqs. (El]). 

It is note that the quark- loop corrections are mode dependent. The assumption of a 
constant gluon invariant mass in FA introduces a large theoretical uncertainty as making 
predictions. In the pQCD approach, however, the gluon invariant mass is related to the 
parton momenta unambiguously and will disappear after the integration. 



D. Magnetic penguins 

This is another kind penguin correction but with the magnetic-penguin operator in- 
sertion. The corresponding weak effective Hamiltonian contains the b -^ sg transition, 

KTf = -%^^^^ts Cl'/Osg, (62) 

with the chromo-magnetic penguin operator, 

Osg = ^m, d, a^^ (1 + 75) ^ G% 6„ (63) 

where i, j being the color indices of quarks. The corresponding effective Wilson coefficient 
Cg = Cgg + C5 [25 . 
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The decay amplitudes obtained by evaluating the Feynman diagrams FigJl^ and FigHb 

can be written as: 

"1 /"OO 

-^3 



'K-qs 



Mi^l = —=Cpm^ / dxidx2dxs / bidbib2db2(j)B{xi,bi) 
v6 Jo Jo 

■{{-{I-X2) [2<Pi{x2) + TK (30^(X2) - <Plix2)) 

+rKX2 (0:^(^2) + 0x(^2))] </>Jfe) 

-r,^ (1 + 0:2) X,<j)i{x2) (3<(X3) + <(X3)) 

-rRTrj, (1 - X2) {(J)k{x2) + 0x(^2)) (30j(x3) - ^J^^l^s)) 

-rKrr,,X3 (1 - 2X2) (0x(^2) - </>^(^2)) (30^(^3) + </)J'^(x3))} 

■Eg{tq)hg{A,B,C,bi,b2,bs,X2) 
- [4:rK^K{x2)(l)^S^3) + 2rxr^,X30K(^2) (30^^ (X3) + (pJ^S^s))] 
■Eg{t'^)hg{A,B\C\b2MM.x^)], 



(64) 



<i. 



:C2m6 



1 ^00 

dxidx2dx^ I bidbib2db2(j)Bixi,bi] 
Jo 



• {{- (1 - X2) [20;^, (X2) + r, (30j;(x2) - 0J,(X2)) 
+r^X2 (0j'^(x2) + 0^,(^2)) (piixs) 

-Tk (1 + X2) X30;J^(X2) (30^(:r3) + 0k(^3)) 

-r,r^ (1 - X2) (0^,(0:2) + 0^^(0:2)) (30^(a;3) - 0^(x3)) 

-r^rxxs (1 - 2x2) (^0^,(^2) - 0^,(^2) j (30;^(x3) + 0x(^3)) 

■Eg(tg)hg{A, B, C, bi, 62, bs, X2) 

- 4r^0^,(^2)0:^(x3) + 2r^rxX30j'^(x2) (30;^(x3) + 0x(^3)) 

■Eg{t'g)hg{A', B', C, 62, 61, 63, 2^1) } . (65) 

Here r^ = ml/mB, r^^ = rriQ/mB. The evolution factors in Eqs. (I64|65p take the form of 



E^^\tj') = alit) elicit) exp[-5„,(t)], 



(66) 



with the Sudakov factor S, 



mg 



Smg{t) = s (ximB/V2, &1 j + s (x2mB/V2, &2J + s n 1 - X2)mB/V2, h 

+S f XsTTls/v^, 63J + S n 1 - X'i)mBl\/2, 63) 



1 

Ti 



In ^^^^/^), + In . ^^(^/^) 



ln(t/A) 
-ln(6iA) ' '"-ln(62A) ' '" - ln(63A) 



In 



(67) 



The hard function hg in the chromo-magnetic penguin amplitude is given by 

p-K/2 

hg{A,B,CMMM.Xi) = -Stixi)KoiBbi)KoiCb3) / rf^tan^ 

Jo 
xJo{Abitan9) Jq {Ab2 tan 9) Jo {Ab^ tan 6*) 



(6^ 
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with the index i = 1,2, the threshold re-summation function St{xi) is given in Eq. ( 1A7I) . 
and 



A = y/x^niB, B = B' = ^/xYX^mB-, C = 2a/(1 - X2) 



XsrriB, 



A' = y^rriB, B' = B, C = ^\xi - x^\mB. (69) 



Here The scale tq, t' , and the gluon invariant mass / and / have been given in Eqs. 
and dSni). 

Finally, the total chromo-magnetic penguin contribution to the considered B — > Kr]^'' 
{K = K^ , K^) decays can be written as 



M^'^ = < K^\H-nB >= -% A, M^^> F,{<P) + M^fJ^ F^ , (70) 



Acmp) T<^r,'\1-/'^'"^P 



x/2 
G 



ria) T^i(A\ , A/rig) 



M™ = < Ki\H-J\B >= -^ A, [M\!l F',{ct>) + Ml^^ F[{ct>)\ . (71) 

The mixing parameters -Fi (</>), F[{(f)), -^2(0) and F2{(j)) have been defined in Eqs. fl22l) and 

ED. 



V. NUMERICAL RESULTS AND DISCUSSIONS 

A. Input parameters 

We use the following input parameters [2|, |40(] in the numerical calculations 

JB = 0.21GeV, Jk = 0.16GeV, m.^, = 547.5MeV, 

rrir,' = 957.8MeV, niK = 0.49GeV, nioK = 1.7GeV, 

Mb = 5.279GeV, m^ = 4.8GeV, Mw = 80.41GeV, 

Tbo = 1.527ps, tb+ = 1.643ps. (72) 

For the CKM quark-mixing matrix, we use the Wolfenstein parametrization as given 
inRef.l,!!^. 

Vud = 0.9745, Vus = \ = 0.2200, IK^I = 4.31 x 10-^, 
Vcd = -0.224, V,d = 0.996, V^ = 0.0413, 
\Vtd\ = 7.4 X 10-^ Vts = -0.042, Vtb = 0.9991, (73) 

with the CKM angles /? = 21.6°, 7 = 60° ± 20° and a = 100° ± 20°. 

B. Branching ratios 

Using the known wave functions and the central values of relevant input parameters, we 
find the numerical values of the corresponding form factors at zero momentum transfer: 

Fi-^{q^ = 0) = 0.37+°:°^(^,), (74) 
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0.40 ± 0.04GeV, which agree well with those obtained in QCD sum rule calcula- 

PP decay can be written as 



for ujb 
tions. 

In the B-rest frame, the branching ratio of a general B 

1 



Br{B -^ M2M3) = Tb 



X \MiB^M2Ms 



(75) 



IGnmB 

where tb is the lifetime of the B meson, x ?« 1 is the phase space factor and equals to 
unit when the masses of final state light mesons are neglected. The total decay amplitude 
in Eq. (175]) is defined as 



M{B^ M2M3) =< M2M3IK// + n 



eff 






(76) 



Using the wave functions and the input parameters as specified in previous sections, 
it is straightforward to calculate the CP-averaged branching ratios for the considered 
four B -^ Ki]^'^ decays, which are listed in Table [B For comparison, we also list the 
corresponding updated experimental results [2] and numerical results evaluated in the 
framework of the QCDF approach 15| . 



TABLE I: The pQCD predictions for the branching ratios (in unit of 10^^). The label LOnlowc 
means the LO results with the NLO Wilson coefficients, and +VC, +QL, +MP, NLO means the 
inclusion of the vertex corrections, the quark loops, the magnetic penguin, and all the considered 
NLO corrections, respectively. 



Mode 


LO 


LOnlowc +VC +QL +MP 


NLO 


Data 


QCDF 


B+ -^ K+r] 


4.7 


4.7 


4.3 4.9 3.1 


3.2 


2.6 ±0.6 


1 q+3.0 
-'^•-'-1.9 


B+ -^ K+n' 


30.2 


46.8 


74.6 48.1 30.2 


51.0 


70.5 ±3.5 


4q 1+45.2 
^^•-'--23. 6 


B^ ^ K% 


3.2 


3.4 


3.1 3.8 2.3 


2.1 


<2.0 


1 1+2-4 


B^ ^ K%' 


3L3 


46.5 


69.7 48.5 20.7 


50.3 


68 ±4 


4fi f^+41.9 



It is worth stressing that the theoretical predictions in the pQCD approach have rel- 
atively large theoretical errors induced by the still large uncertainties of many input 
parameters, such as quark masses(m„^rf,ms), chiral scales (mox, wig, mp), Gegenbauer co- 
efficients (Oj. '''%■■■), LUh and the CKM angles (0,7), etc. The NLO pQCD predictions 
for the CP-averaged branching ratios with the major theoretical errors are the following 



Br{ B+ ■ 

Br{ 5+ - 

Br{ 5° 

Br( B° - 



K+T]') 

^ir%) 



+1.2 



13.5 



[3.2 

[51.0_8.2 

I- 
[50.3+^^2^ 



il?i( 



m. 



);}:^( 



i^)] X 10- 



UJh 



+11.2, 
'-6.2 I 



m. 



[2.lt'oii^,)tliim, 



1+1.0/ »? 
'-0.9l"2 



X 10" 



X 10" 



Hni'2'^ 



m. 



)-2.7('^2 



X 10" 



0.4±0.04GeV, m. 



(77) 
(78) 
(79) 
(80) 

130 ±30 



The major errors are induced by the uncertainties of ujb 

MeV and Gegenbauer coefficient al = 0.44±0.22 (here a^ denotes a^'' or 0^2), respectively. 

In Figs. [8] and [H we show the parameter dependence of the pQCD predictions for the 
branching ratios of 5+ -^ K^r]^'^ and B^ — > K^r]^'^ decays for cj^ = 0.4 ± 0.04 GeV, 
7= [0°,180°]. 

From the numerical results about the branching ratios, one can see that 
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FIG. 8: The 7 dependence of the branching ratios (in units of 10 ^) of B^ -^ K^rj^'> decays for 
cob = 0.36 GeV (dotted curve), 0.40 GeV (sohd curve) and 0.44 GeV (dashed curve). 



The decay amplitude B -^ Krjq and B —>■ Krjs interfere constructively for B -^ Kr]' 
decays, but destructively for B —>■ Krj decays. This mechanism results in a factor 
of 6 — 10 disparity for the branching ratios of -B — * K^rj' and B — > K^rj decays. 

The LO pQCD predictions for branching ratios are much smaller (larger ) than the 
measured values for B -^ Kr]' [B —>■ Krj) decays, show the same tendency as found 
in Ref. fTT 



The NLO contributions can interfere constructively (destructively) with the cor- 
responding LO parst for B -^ Krj' ( B —>■ Krj) decays. For B^ —>■ K^rj' and 
B'^ -^ K^T]' decays, the NLO contributions provide a 70% enhancement to their 
branching ratios . For B^ —>■ K^rj and B~^ — > K~^ri decays, on the other hand, the 
NLO contributions give rise to a 30% reduction to their branching ratios and result 
in the good agreement between the pQCD predictions and the data. 

The NLO pQCD predictions for branching ratios Br(B — > Krf^'^) agree very well 
with the measured values within one standard deviation. The NLO contributions 
play an important role in understanding the observed pattern of branching ratios 
of the four B -^ Kr]^'^ decays. 



C. CP- violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries oi B ^ Kr]^'^ decays 
in pQCD approach. For B^ -^ K^r]^'' decays, the direct CP-violating asymmetries Acp 
can be defined as: 



Adir \_ \ 



Mf\'-\Mf\' 



^CP 



\Mj 



\^ + \Mf\^' 



(81) 
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FIG. 9: The 7 dependence of the branching ratios (in units of 10 ^) of B^ — > K'^r]^'' decays for 
ujb = 0.36 GeV (dotted curve), 0.40 GeV (sohd curve) and 0.44 GeV (dashed curve). 



Using Eq. (IHTj) . it is easy to calculate the direct CP-violating asymmetries for the 
considered decays, which are listed in Table [TTl As a comparison, we also list currently 
available data [^ and the corresponding QCDF predictions [l5|]. 



TABLE II: The pQCD predictions for the direct CP asymmetries in the NDR scheme (in units 
of 10^^), the QCDF predictions 15|] and the world average as given by HFAG [4]. 



Mode 


LO 


^Onlowc +VC +QL +MP 


NLO 


Data 


QCDF 




9.3 
-10.1 


10.3 31.1 7.8 7.6 
-7.3 -10.6 -5.9 -10.4 


-11.7 
-6.2 


-27 ±9 
1.6 ±1.9 


-,0 Q+29.0 



The NLO pQCD predictions for A'^Spi.B^ 
theoretical errors are 



K^r]^'^) (in unit of 10 ^) with the major 



Atp{B^ 
Atp{B^ - 






-ii-7+^:6K 



^+2-9(a:?^^ 



-4.2I /y'-5.6V"'2 

-6.2ll:?(m.)li:g(7);}:^(a^'), 



(821 



where the dominant errors come from the variations oirris = 130 ±30 MeV, 7 = 60° ±20° 
and Gegenbauer coefficient a^'' = 0.115 ± 0.115, respectively. 

As to the CP-violating asymmetries for the neutral decays B^ -^ K^r]^'\ the effects of 
B^ — B^ mixing should be considered. The CP-violating asymmetry of B^(B^) -^ K'^r]^'^ 
decays are time dependent and can be defined as 



P i?o(At) ^fcp)-T (i?0(At) ^ fcp) 



Acp 

P (i?o(At) -. fcp) + V (50(At) ^ fcp) 

= A'^jp cos(AmAt) ± A^]? sin(AmAt), (83) 

where Am is the mass difference between the two B^ mass eigenstates. At = tcp — Uag 
is the time difference between the tagged B^ {B ) and the accompanying B {B^) with 
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opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct and 
mixing induced CP-violating asymmetries A^p ( or Af in term of Belle Collaboration) 

mix 
CP 



and ATi'^S can be written as 



A' 



dir 
CP 



A 



lA 



CP\ 



1 



/ 



l + |Acp| 
with the CP-violating parameter \cp 



A mix 
•^CP 



Sf 



2/m(A 



CP) 



l + IA 



CP\ 



U) 



A, 



CP 



{fcp\Heff\B ) 
{fcp\H,ff\BO)' 



By integrating the time variable t, one finds the total CP asymmetries for B^ 

decays, 



^5) 



K%(') 



A tot 

■/\cp 



l + x- 



A dir I 
.^CP + 



X 



l + x 



A mix 
2^CP ' 



(86) 



where x = Am/T = 0.775 [4C 

In Table [nil we show the pQCD predictions for the central values of the direct, mixing- 
induced and total CP asymmetries for B^ — >■ K^t]^'^ decays, obtained by using the LO or 
NLO Wilson coefficients, and adding the vertex corrections, the quark loops, the magnetic 
penguin, or include all the mentioned NLO corrections, respectively. 



TABLE III: The pQCD predictions for the direct, mixing induced and total CP asymmetries 
(in units of 10^^) for B^ -^ K^rj^'' decays, and the world average as given by HFAG [^. 



Mode 


LO 


^Onlowc +VC +QL +MP 


NLO 


Data 


Afp{B^ - K^si) 


-4.2 

1.4 


-1.5 -11.2 -0.9 -1.9 
0.0 1.5 0.7 -0.1 


-12.7 
2.3 


9±6 




61.6 

64.6 


67.3 64.4 66.9 67.9 
63.5 63.4 63.2 63.2 


61.9 

62.7 


61 ±7 


A'SUB'^Klr^) 


27.2 
32.1 


31.7 24.2 31.9 31.7 

30.8 31.7 31.0 30.5 


22.1 
31.8 






The NLO pQCD predictions for ^gj,(5° - 
of 10~^) with the major theoretical errors are 



.^0 

A^piB 



A dir 

•^CP( 

\dir 







Amix ( dO 

cp{B 



^CP{B 

Amix f jjO 






/sfOr^W) and A'S'piB^ 



-12.7±4.1(m.)+?;^(7)^f,(a^^), 



2.3^°i(m.)+°:^(7)^°:?(a^'), 



Q^ g+35.8,'„,^+35.3, 



-65.0 



(7): 



■64.3 



62.7tllih)^''-'' 



-64.7 



a) 



a) 



K^7]^"i) (in unit 



^7) 



where the dominant errors come from the variations of m^ = 130 ±30 MeV, 7 = 60° ±20°, 
a = 100° ± 20°, and the Gegenbauer coefficient a^'' = 0.115 ± 0.115, respectively. 
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FIG. 10: The NLO pQCD predictions for direct CP asymmetries (in percentage) oi B^ -^ K^r]^'' 
and B — > K r]^'' decays. 
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FIG. 11: The 7-dependence (a) and the a-dependence (b) of the total CP-asymmetries of 
B^ -^ KgT] (sohd curve) and B^ -^ K^rj' (dotted curve) decays. 



In Fig.lini we shown the 7-dependence of the pQCD predictions for direct CP- violating 
asymmetries of 5° — > K'^t]^'^ and B~^ — >■ K~^ri^'^ decays. In Fig. [TT|, we shown the a- 
dependence of the total CP-violating asymmetries for B^ -^ Kgi] (solid curve) and B^ -^ 
KgT]' (dotted curve), respectively. 

From the pQCD predictions and currently available experimental measurements for 
the CP violating asymmetries of the four B -^ Kr]^'^ decays, one can see the following 
points: 

(a) For B^ — » K^t] decay, the measured direct CP asymmetry is 3 standard devi- 
ation from zero. The LO pQCD prediction changed its sign and become consistent 
with the measured one due to the inclusion of NLO contributions. 

(b) For A'^'piB^ -^ K^V'), the pQCD prediction is changed from -10% to -6% 
due to the inclusion of NLO contributions, which is consistent with the measured 
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zero result within one standard deviation. 

(c) For B^ -^ K^T]^'^ decay, the effects of NLO contributions to their CP asymmetries 
are rather small, as can be seen from the numerical results as given in Table IIIII 

(d) For neutral fi° -^ K^t]^'^ decays, the PQCD predictions are A'^'p{B° -^ K^t]') ^ 
2.3% and A'SpiB^ -^ K^sV') ~ 63%, which agree very well with the data: (9 ± 6)% 
and (61 ± 7)%. This means that the deviation AS = —TjfSf — sin 2(3 for B^ —>■ K^rj' 
decay is also very small in the pQCD approach. 



VI. SUMMARY 

In this paper, we calculated the branching ratios and CP-violating asymmetries of 
B^ -^ K~^r]^'^ and 5° -^ K^r]^'' decays in the pQCD approach. The partial NLO con- 
tributions considered here include: QCD vertex corrections, the quark-loops and the 
chromo-magnetic penguins. 

From our calculations and phenomenological analysis, we found the following results: 

(a) The pQCD predictions for the form factors of -B ^ r/'^'-' and B ^ K transitions 



are 



F^^'\q^ = 0) = 0.21 ± 0.03(tUb), 
F^~*'^\q^ = Q) = 0.17±0.03(cub), 



for ujb = 0.40 ± 0.04 GeV, which agree well with those obtained in QCD sum rule. 

(b) For branching ratios, the NLO pQCD predictions (in unit of 10^^) are 

Br{B+ -^ K+r]) = 3.2+1% 

Br{B+ - K^) = 51.0™, 

5r(5° ^ K'^T]) = 2.ltli, 

Br{B^ -^ K^^r]') = 50.3+}[;:^, (89) 

where the individual theoretical errors have been added in quadrature. The decay 
amplitude B -^ Krjg and B -^ Krjg interfere constructively for B -^ Kt]' decays, but 
destructively for B —>■ Krj decays. The NLO contributions in the pQCD approach, 
furthermore, can provide a 70% enhancement to Br{B — > Krj'), but a 30% reduction 
to Br{B — > Kt]). The large branching ratio oi B —>■ Krj' decays, as well as the large 
disparity Br{B -^ Kt]') ^ Br{B -^ Kif) can therefore be understood naturally. 

(c) The pQCD predictions for the CP asymmetries of -B ^ Krj^''' decays are con- 
sistent with currently available data. For neutral B'^ — > K'^r]^'^ decays, for example, 
the PQCD predictions are ^^'^(5° -^ K^s^l') ~ 2.3% and A'g^iB^ -^ K'^gVl ~ 63%, 
which agree very well with the measured values of (9 ± 6)% and (61 ± 7)%, respec- 
tively. 
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(d) In this paper, only the partial NLO contributions have been taken into account. 
We think that they are the dominant part of the whole NLO corrections. To achieve 
a complete NLO calculations in the pQCD approach, the still missing pieces from the 
emission diagrams, hard-spectator and annihilation diagrams, should be evaluated 
as soon as possible. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function /ij's, coming from the Fourier transformations of the hard 
kernel H^^\xi, bi), 

he{xi,X2, h, 62) = Ko (y/x^niBbi) [9{bi - b2)Ko {^/x^niBbi) Iq {^/x^mBb2) 

+0{b2 - bi)Ko iy^mBb2) h (v^"^b^i)] ^'((xs), (Al) 



ha{x2, X3, 62, &3) = Ko \i\/{l - X2)x^mBb2] [9{b^ - &2)-f^o (iy/xsniBb^) Iq {iyJx^mBb2) 

+0{b2 - 63)^0 (^v^"^s^2) ^0 (iy^ruBbs)] Stix^), (A2) 



hf{xi,X2,Xs,bi,b3) = <9{bl-b3)Ko{MB^/x^bl)lo{MB^/x^b 
+ e{bs - 6i)Ko(Mbv/SI^63)Io(Mbv/SI^6i) 



fUoi^/ix^ix^^^x^MBbs), for xi - x^ < 
Kq (V(a;2(a;i - X3)MBb3), for xi - x^ > 



hl{xi,X2,X3,bi,b3) = \ 9{bi - 63)Ko(iA/(l - X2)x3biMB)lo{i\/{l - X2)x3b3ME 



(A3) 



B, 



+ {9{bs - 6i)Ko(V(l - X2)x3b3MBMtVi^ - X2)x3blM, 

Ko{MB^/(x^'^^X3)(T^^X2jbl), for xi - ^3 > . .^^. 



2 



Hq \MB^/{x^^^Xl}(J^^^X2}bl), for X1-X3 < 
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hj{xi,X2,X3,bi,b2) = < 9{hi - 63)Ko(i V(l - X2)x'ihiMB)lo{i\/ {I - X2)xsb3MB) 

Ko(MbFi6i), for F^ > o' 

r(i 



f}l','>{MB./\Ff\bi), for F2<0^ 



(A5) 



where Jq is the Bessel function and Kq, Iq are modified Bessel functions with KQ{—ix) 
— {n/2)Yo{x) + i{TT/2)Jo{x), and i^(i)'s are defined by 



"h 



1 -X2(l - Xi -X3). 



(A6) 



The threshold resummation form factor St{xi) is adopted from Ref. [4l|. It has been 
parameterized as 






(AT) 



where the parameter c = 0.3. This function is normahzed to unity. 

The evolution factors Eg , and Ea , appeared in the decay amplitudes are given by 

Ee{t) = as(t)exp[-Sab{t)], 

E'^{t) = as{t)exp[-Scd{t)]\b2=bi, 

Kit) = as{t)exp[-Sef{t)]\b2=b3, 

Ea{t) = asit)exp[-Sghit)]. {Ai 

The Sudakov factors used in the text are defined as 

Sab{t) = s (ximB/V2, bij + s (x2mB/V2, 62J + s N 1 - X2)mB/V2, h 



1 

"a 



In 



ln(t/A) 



In 



ln(t/A) 



-ln(6iA) -ln(62A) 

Scd{t) = s (ximB/V2, hi) + s (x2mB/V2, hj + s Ul - X2)mB/V2, h 

+s (x3mB/V2, 63J + s n 1 - X3)mB/V2, k 
1 



(A9) 



Pi 



21nJ^^(^ + ln.l^(^/^) 



ln(6iA) -ln(63A) 

Sefit) = s \^XimB/V2,bi] + s (x2TnB/V2,h] + s Ul- X2)mB/V2,b'- 

+s \^X3mB/V2, 63J + s (^(1 - X3)mB/V2, 63J 

_J_r ln(t/A) ln(t/A) 

A [ -ln(6iA)^ -ln(62A) 

Sgh{t) = s (x2mB/V2, 62] + s [x^niBlV^, b 
+ s\{l- x-i)mBl\f2, 63 j - — 



(AlO) 



s nl -X2)mBlV2,br^ 



l^J^(VAl_^l^,ln(t/A) 



ln(63A) 



ln(62A) 



(All) 



(A12) 
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where the function s{q, b) are defined in the Appendix A of Ref. 29|. The scale t^'s in the 
above equations are chosen as 



t'a 

h 
A 

tc 

td 



max 
max 
max 
max 
max 
max 
max 
max 



X2mB, ^/x^ruB, 1/&1, 1/&2) , 

xTms, y/x^niB, l/bi, I/62) , 
a/x2|1 -X3 -Xi\mB, ^/x^TJiB, l/bi, 1/63) , 
V^X2|x3 -Xi\mB, y/x^niB, 1/61, 1/&3) , 
a/(1 -X2)x3mB, a/|xi -X3|(1 -X2)mB, l/bi, 1/63) , 
a/|1 -X2(1 - X3-Xi)\mB, a/(1 - X2)x3mB, 1/61, 1/63) 
a/(1 -X2)x3mB, a/(1 -X2)mB, 1/&2, I/&3) , 
V/(1 - X2)x2.mB, y/xiniB, 1/&2, I/&3) • 



(A13) 
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